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The random field Ising model is studied numerically at both zero and positive temperature. 
Ground states are mapped out in a region of random and external field strength. Thermal states and 
thermodynamic properties are obtained for all temperatures using the the Wang-Landau algorithm. 
The specific heat and susceptibility typically display sharp peaks in the critical region for large 
systems and strong disorder. These sharp peaks result from large domains flipping. For a given 
realization of disorder, ground states and thermal states near the critical line are found to be strongly 
correlated-a concrete manifestation of the zero temperature fixed point scenario. 
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The random field Ising model (RFIM) is one of the sim- 
plest non-trivial spin models with quenched disorder. De- 
spite thirty years of study it is still not well understood. 
It has been proved that an ordered phase exists for suffi- 
ciently low temperature and dimension d > 2 0, 0, 0, ■ 
The phase transition between the ordered and disordered 
phases for d > 2 is believed to be continuous and con- 
trolled by a zero temperature fixed point 0, El 0} ■ Cur- 
rently, there is no controlled renormalization group anal- 
ysis of the RFIM phase transition and Monte Carlo sim- 
ulations [il 0, 0, 0] are restricted to small systems and 
have been inconclusive. As the strength of the random 
field increases the transition moves to lower temperature 
and the critical line intersects the zero temperature line 
at a zero temperature phase transition. Numerical stud- 
ies of the zero temperature transition 0,0,0,0 play 
an important role in understanding the model. Ground 
states are much easier to simulate than thermal states 
and, according to the zero temperature fixed point hy- 
pothesis, the T = and T > transitions are in the same 
universality class. Critical exponents have been obtained 
from zero temperature studies that are mostly consistent 
with the scaling theories 0, 0] , series methods and 
real space renormalization group approaches [TtI list Il9l | . 

In this paper wc present numerical results at both 
T = and T > for the same realizations of ran- 
dom fields. For T > we use the Wang-Landau 
and Metropolis algorithms. For T = we find ground 
states using the push-relabel algorithm 00 A major 
conclusion of the paper is that spin configurations found 
near the critical line are strongly correlated with ground 
states near the zero temperature critical point. This ob- 
servation is consistent with the original Imry-Ma anal- 
ysis, incorporated in the zero temperature fixed point 
scenario, that the large scale properties of the critical 
point depend on the competition between random fields 
and couplings with thermal fluctuations serving only to 
renormalize the strength of these couplings. However, 
the correlation found here for single realizations of disor- 
der along the critical line is not implied by the existence 
of a zero temperature fixed point, which implies only the 



similarity of zero temperature and positive temperature 
critical ensembles. 

The Hamiltonian of the RFIM studied in this paper is 



n 
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The summation is over all nearest neighbors i and 
j on a simple cubic lattice with periodic boundary con- 
ditions, spins Si take the value ±1, A is the strength of 
disorder, hi is the random field chosen from a Gaussian 
distribution with mean zero and variance one, and H is 
the external field. Two important quantities are the mag- 
netization (order parameter) m — (1/ L 3 )^ li and the 
bond energy e = (1/L 3 ) j\ SiSj. We define the dis- 
order strength separately from the normalized random 
fields because one of our primary concerns is to examine 
single realizations of random fields as disorder strength, 
temperature and external field are varied. Previous an- 
alytic [22| and real space renormalization group stud- 
ies ^3 also considered single realizations of disorder at 
the phase transition but do not compare realizations at 
different disorder strengths as is done here. 

Consider the set of ground states of a single realization 
of disorder. We obtain these using a method first intro- 
duced by Ogielski 0] . To determine the ground state at 
a given value of H and A, the RFIM problem is mapped 
onto the MAXFLOW problem, which is then solved us- 
ing the push-relabel algorithm 000 0. The set 
of all ground states in a region in the H — A plane is 
mapped out using a method described in [24| and simi- 
lar to the techniques discussed in 0, 0] . Figure ^ is 
a portrait of all the ground states of a single realization 
of random fields in a 32 3 system in a small region of the 
H — A plane near the finite size critical point, discussed 
below. Each line represents values of the parameters for 
which two ground states are degenerate and across each 
line a single connected domain is flipped. Within each 
polygon bounded by these lines, a single spin configura- 
tion is the ground state. At points where two lines cross, 
four ground states are degenerate and the four configu- 
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rations differ by the orientation of two separate domains. 
More interesting are "triple points" where a line bifur- 
cates into two lines in a Y shape. At triple points three 
ground states are degenerate but the three domains cor- 
responding to the three lines are not independent. The 
spin configuration at the top of the Y results from the 
break-up of the large domain that flips across the verti- 
cal line of the Y as shown schematically in Fig. [3 The 
triple point has some characteristics of a thermal first- 
order transition where two ordered states co-exist with a 
disordered state. 

When a coexistence line is crossed and a domain is 
flipped physical quantities except for the total energy are 
discontinuous. To visualize the size of the discontinuity, 
lines are drawn with a thickness that is proportional to 
the jump in the magnetization. The picture is simplified 
by removing the large number of lines with small bond 
energy jump (Se < 0.3), as shown in Fig.^>. The simpli- 
fied picture reveals a tree-like structure built from triple 
points. The triple point with the largest energy disconti- 
nuity is located at the center of the picture. In the region 
above this triple point the magnetization is small while 
the line extending below the triple point is the coexis- 
tence line separating the plus and minus ordered states. 
In Ref. this triple point was identified as the finite- 
size critical point and its scaling properties were studied. 
The size of the discontinuity in the bond energy is gov- 
erned by the specific heat exponent. We have also exam- 
ined the large discontinuities in bond energy and magne- 
tization along the H = and shown j3| that these scale 
with the specific heat exponent and magnetic exponents, 
respectively. Within a region that shrinks as L^ a+I3 ^ 2 ^ 1 ' 
and L x l v in the H and A directions, respectively, the 
tree-like structure is statistically self-similar but not self- 
averaging-each realization has a unique tree-like struc- 
ture. 

We study the RFIM as a function of temperature us- 
ing the Wang-Landau [lOl and the Metropolis algorithms. 
The Wang-Landau algorithm is a flat histogram Monte 
Carlo method that automatically determines the density 
of states. Thermodynamic quantities at all temperatures 
are then derived from the density of states and the statis- 
tics of the magnetization as a function of energy. The al- 
gorithm smooths the energy landscape and is much more 
efficient than the conventional Metropolis algorithm for 
sweeping a range of temperatures. Once a temperature 
is chosen for detailed study, the Metropolis algorithm is 
used to find the thermally averaged spin configuration. 
We determined the specific heat and susceptibility for 
systems up to size 32 3 . We find that for large enough 
systems (> 16 3 ) and strong enough disorder, the spe- 
cific heat and the susceptibility typically display one or 
more sharp peaks. In Fig. |21 we show the specific heat 
and the susceptibility as a function of temperature for 
the same realization of normalized random fields whose 
ground states are shown in Fig. ^ The random field 



strength is Ao = 2.0 and the external field is set to zero. 
Two sharp peaks appear in both quantities at the same 
temperatures. We have simulated one hundred 16 3 real- 
izations with A = 1.5 and find that about 1/3 of them 
have sharp peaks. The number increases to 1/2 if the 
random field is strengthened to Ao = 2.0. For size 32 3 
and Ao = 2.0 we have simulated nine realization and 
sharp peaks are observed for all of them. We tentatively 
conclude that the probability of sharp peaks appearing 
increases with the system size and the strength of random 
field. 

The sharp peaks in the specific heat and susceptibil- 
ity can be understood within the zero temperature fixed 
point picture of the RFIM phase transition. This pic- 
ture predicts that the behavior in the critical region at 
finite temperature is determined by the competition be- 
tween couplings and random fields with thermal fluctu- 
ations serving only to renormalize the strength of these 
quantities. One conclusion of this paper is that this sce- 
nario appears to be true for individual realizations of 
normalized random fields. The sharp peaks in the ther- 
modynamic quantities can be matched one to one with 
the large jumps at zero temperature. Furthermore, the 
spin configurations on either side of the sharp peaks can 
be mapped onto the ground states on either side of the 
corresponding large jumps. 

For a single realization of random fields, we obtain the 
thermally averaged spin configuration near the peaks at 
finite temperature, and compare these thermal states to 
the ground states near the two largest jumps at zero tem- 
perature. Figures^Ji, e and f show one plane through the 
system with Ao = 2.0 and at temperatures just before 
peak 1 (T = 2.2), just after peak 1 (T = 2.5), and just af- 
ter peak 2 (T = 2.8), respectively. The difference among 
the states shows that the sharp peak corresponds to flip- 
ping a relatively large domain. It is evident that these 
three states are strongly correlated with the ground state 
spin configuration before the jump 1 (A = 2.36), just af- 
ter jump 1 (A = 2.41), and just after jump 2 (A = 2.54), 
as shown in Fig. 0^,, b and c, respectively. (The labels of 
jumps and peaks are given in Fig.^> and Fig. El) Simi- 
lar correlations between ground states and thermal states 
were found in one dimension 

El- 
Some correlation between ground states and thermal 
states persists to much smaller values of Ao in a regime 
where the thermodynamic properties no longer display 
sharp peaks. Figure^, h and i show the same realization 
of disorder and the same plane through the system but 
with Ao = 0.5. Here the specific heat has a rounded 
peak at T — 4.375. Figures 0k, h and i correspond to 
temperatures 4.0, 4.3 and 4.45, respectively. Although 
there is considerable thermal "blurring" in these pictures, 
evidence of the ground state is unmistakable. 

A quantitative characterization of the correlation be- 
tween ground states and thermal states for the same re- 
alization can be obtained from the correlation measure, 
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FIG. 1: Ground states of the RFIM in the H — A plane, (a) All the ground states of a single 32 3 realization of disorder. Along 
each line two ground states coexist that differ by flipping a single connected domain. The thickness of a line is proportional to 
the magnetization jump across the line.(b) The same realization as in (a), but only lines with the bond energy jump Se > 0.03 
are shown. Along the H = axis there are two major jumps, which are labeled as 1 and 2 in the graph. 
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FIG. 2: Schematic picture of a triple point. The shaded ovals 
show the orientation of spins within a single domain that flips 
crossing the vertical line and is broken into two pieces crossing 
the diagonal lines. 
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where the overbar is an average over realizations of disor- 
der and (sJa.t is the thermal average of the spin at the 
ith site at disorder A and temperature T or, if T = 0, 
it is the ground state spin value. For each realization, 
the temperature T* = T max + 0.1 where T max is the tem- 
perature of the maximum of the specific heat; one of the 
sharp peaks in C if sharp peaks exist. Thus, for each 
realization, we pick a thermal state just above the tran- 
sition temperature. Figure [S] shows q vs. A for sizes 16 3 
and 32 3 and A = 1.5, with 96 realizations for size 16 3 
and 9 for size 32 3 . A peak in the correlation occurs at 
A w 2.65 where q rs 0.75. The value, A « 2.65, is about 
0.15 larger than the average A at the largest disconti- 
nuity in the bond energy for system size 32 3 . The inset 
in Fig. shows the average correlation between thermal 
states of one realization and ground states of another 
for size 16 3 , which is nearly zero as expected. A second 
measure, q* is obtained by choosing the value A* for each 
ground state realization to give the largest correlation to 
the thermal state and then averaging over realizations. 
We find that for size 32 3 , q* = 0.80 ± 0.06 for A = 1.5 
and q* = 0.85 ±0.05 for A = 2.0. Together, these result 
provide quantitative confirmation that the thermal states 
at temperatures slightly above the thermal critical point 



are strongly correlated with the ground states at disor- 
der strength slightly higher than the zero temperature 
critical point. 

The strong correlations between states at different 
temperatures is ostensibly in conflict with the idea of 
"chaos" in the RFIM. Chaos in systems with quenched 
disorder, such as spin glasses and the RFIM, refers to the 
sensitivity of spin configurations to small perturbations 
either in temperature or in quenched disorder [lH Erj Ejf . 
The existence of chaos in the RFIM is controversial and 
is not definitively established. This work suggests that 
chaos is not present along trajectories in the A — T plane 
following the critical line. 

In summary, we find that sharp peaks in thermody- 
namic functions resulting from the flipping large do- 
mains are typical near the critical point. In addition, 
spin configurations near the transition are similar to the 
ground states near some corresponding large jump at zero 
temperature. If this connection between critical ground 
states and thermal states persists to large system size it 
supports a strong version of the zero temperature fixed 
point scenario: the sequence of states near the zero tem- 
perature critical point obtained by varying A for T = 
can be mapped onto the sequence of thermal states near 
the critical point obtained by varying T for fixed values 
of A , A < A c . 

We thank Alan Middleton for useful discussions. This 
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FIG. 3: The specific heat and the susceptibility of the same 32 3 realization as in Fig. Qwith A = 2.0 and H = 0. Two sharp 
peaks, labeled 1 and 2, are observed, which correspond to the two large jumps 1 and 2 in Fig.0 respectively. 




FIG. 4: Spin configurations near the critical points at zero temperature and finite temperatures for a single realization of 
normalized random fields. Each panel is the same plane of a 32 3 realization with black representing spin down; white, spin 
up; and shades of gray, the thermally averaged spin state. From left to right in the top two rows, panels are at A (T) before, 
between and after jumps (peaks) 1 and 2 in Fig. (Fig. |3J. Specifically, panels (a), (b) and (c) are ground states at A = 2.36, 
2.41 and 2.54, respectively. Panels (d), (e) and (f) are at A = 2.0 and T — 2.2, 2.5 and 2.8, respectively. Panels (g), (h) and 
(i) are at A = 0.5 and temperatures 4.0, 4.3 and 4.45, near the peak in the specific heat at T = 4.375. 
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FIG. 5: Disorder averaged correlation q of a thermal state just above the transition temperature at Ao = 1.5 to ground states 
at disorder stength A for the same realization of random fields. Solid squares for size 16 3 and open circles for size 32 3 . Only a 
few error bars are drawn to make the figure easy to read. The inset shows the correlation of thermal states with ground states 
of a different random field realization. 
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